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Prologue 

Background. The hehcity of a magnetic field or fluid flow measures the extent to which its 
orbits wrap and coil around one another. Woltjer introduced this notion in the late 1950s during 
his study of the magnetic field in the Crab Nebula, showed that the helicity remains constant as the 
field evolves according to the equations of ideal magnetohydrodynamics, derived from this a lower 
bound for the changing field energy, and calculated the stable field at the end of the evolution. 
The term "helicity" was coined ten years later by Moffatt, who rewrote Woltjer's integral formula 
to reveal its analogy with Gauss's linking integral for two disjoint closed curves in 3-space. 

When the helicity of a magnetic field or fiuid fiow is zero, the lower bound for energy that it 
provides is also zero, and one hopes for a "higher order helicity" which can provide its own lower 
bound for energy. Monastyrsky and Retakh [1986] and Berger [1990] prepared the way for this via 
integral formulas derived from the Massey product formulation of Milnor's triple linking number of 
a three-component link, and various authors since then have shown in special cases how to derive 
nonzero lower energy bounds. 

What we do here. We describe a new approach to triple linking invariants and integrals, 
aiming for a simpler, wider and more natural applicability to the search for higher order helicities. 

To each three-component link in Euclidean 3-space, we associate a generalized Gauss map 
from the 3-torus to the 2-sphere, and show that the pairwise linking numbers and Milnor triple 
linking number that classify the link up to link homotopy correspond to the Pontryagin invariants 
that classify its generalized Gauss map up to homotopy. 

When the pairwise linking numbers are all zero, we give an integral formula for the triple 
linking number analogous to the Gauss integral for the pairwise linking numbers, but patterned 
after J.H.C. Whitehead's integral formula for the Hopf invariant, and hence interpretable as the 
ordinary helicity of a related vector field on the 3-torus. 

What's new about this? Our generalized Gauss map from the 3-torus to the 2-sphere is 
a natural extension of Gauss's original map from the 2-torus to the 2-sphere; like its predecessor 
it is equivariant with respect to orientation-preserving isometrics of the ambient space, attesting 
to its naturality and positioning it for application to physical situations; it applies to all three- 
component links, not just those with pairwise linking numbers zero; and when the pairwise linking 
numbers are zero, it provides a simple and direct integral formula for the triple linking number 
which is a natural successor to the classical Gauss integral, with an integrand invariant under 
orientation-preserving isometrics of the ambient space. 

Application. Komendarczyk [2009, 2010[ has applied this approach in special cases to derive 
a higher order helicity for magnetic fields whose ordinary helicity is zero, and to obtain from this 
nonzero lower bounds for the field energy. 

In the first paper of this series [2011[, hereafter "Part 1", we did all of the above for three- 
component links in the three-sphere. The first step there was to find a geometrically natural 
generalized Gauss map. That same first step is taken here in Euclidean 3-space, but the map itself 
is entirely different because the requirement of geometric naturality involves a different, and in this 
case non-compact, group of isometrics. After describing this new version of the generalized Gauss 
map, we build a bridge between the spherical and Euclidean versions, across which we transport 
proofs and save labor. 



1 



DeTurck, Gluck, Komendarczyk, Melvin, Nuchi, Shonkwiler and Vela-Vick 



I. Introduction 

Setting the stage. 

Three-component links in Euclidean 3-space were classified up to link homotopy - a de- 
formation during which each component may cross itself but distinct components must remain 
disjoint - by John Milnor in his senior thesis, published in 1954. A complete set of invariants is 
given by the pairwise linking numbers p , q and r of the components, and by the triple linking 
number, which is the residue class /i of one further integer modulo the greatest common divisor 
of p, q and r. 

For example, the Borromean rings shown below have p — q — r — Q and /x — ±1, where the sign 
depends on the ordering and orientation of the components. 




Borromean Rings 

This is a photograph, courtesy of Peter Cromwell, of a panel in the carved walnut doors of the 
Church of San Sigismondo in Cremona, Italy. 

To each ordered, oriented three-component link L in M-^, we will associate a generalized Gauss 
map gL from the 3-torus = x x to the 2-sphere , in such a way that link homotopies 
of L become homotopies of g^. The definition of will be given below. 

Maps from to were classified up to homotopy by Lev Pontryagin in 1941. A complete set 
of invariants is given by the degrees p, q and r of the restrictions to the 2-dimensional coordinate 
subtori, and by the residue class v of one further integer modulo twice the greatest common divisor 
of p , q and r, the Pontryagin invariant of the map. 

This invariant is an analogue of the Hopf invariant for maps from S'^ to , and is an absolute 
version of the relative invariant originally defined by Pontryagin for pairs of maps from a 3-complex 
to the 2-sphere that agree on the 2-skeleton of the domain. 

Our first main result. Theorem A below, equates Milnor's and Pontryagin's invariants p, q and 
r for L and gL, and asserts that 

2fi{L) = 

In the special case when p = g = r = 0, we derive an explicit and geometrically natural integral 
formula for the triple linking number, generalizing Gauss's classical integral formula for the pairwise 
linking number and patterned after J.H.C. Whitehead's integral formula for the Hopf invariant. 
This formula and variations of it are presented in Theorem B below. 

In the rest of this introduction, we give the background and motivation for our work, then lead 
up to and provide the definition of the generalized Gauss map of a three-component link in M.^, give 
careful statements of Theorems A and B, and finally present the results of a numerical calculation 
of Milnor's triple linking number using Theorem B. 
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Background and motivation. 

We recall the famous integral formula of Gauss [1833] for the linking number of two disjoint 
smooth closed curves 

X = {x{s) : seS^} and Y = {y{t) : t e S^} 

in Euclidean 3-space M'^: 

^ 4tt Jj,2 ds dt \x - 

The helicity of a vector field V defined on a bounded domain il in is given by the formula 

Hel(y) = / V{x) X V{y) • dx dy , 

where dx and dy are volume elements. 

There is no mistaking the analogy with Gauss's linking integral, and no surprise that helicity is 
a measure of the extent to which the orbits of V wrap and coil around one another. 

Woltjer [1958] introduced this notion during his study of the magnetic field in the Crab Nebula, 
showed that the helicity of a magnetic field remains constant as the field evolves according to the 
equations of ideal magneto-hydrodynamics, derived from this a lower bound for the changing field 
energy, and calculated the stable field at the end of the evolution. The term "helicity" was coined 
by Moffatt ]1969], who also derived the above formula from Woltjer's original expression. 

Since its introduction, helicity has played an important role in astrophysics and solar physics, 
and in plasma physics here on earth. 

Our study was motivated by a problem proposed by Arnol'd and Khesin [1998] regarding the 
search for "higher helicities" for divergence-free vector fields. In their own words: 

The dream is to define such a hierarchy of invariants for generic vector 
fields such that, whereas all the invariants of order < k have zero value 
for a given field and there exists a nonzero invariant of order A; + 1, this 
nonzero invariant provides a lower bound for the field energy. 

Previous integral formulas for Milnor's triple linking number and attempts to define a higher 
order helicity can be found in the work of Massey [1958, 1969], Monastyrsky and Retakh ]1986], 
Berger ]1990, 1991], Guadagnini, Martelfini and Mintchev [1990], Evans andBerger [1992], Akhmetiev 
and Ruzmaiken [1994, 1995], Arnol'd and Khesin [1998], Laurence and StreduHnsky [2000], Leal [2002], 
Hornig and Mayer [2002], Riviere [2002], Khesin [2003], Bodecker and Hornig [2004], Auckly and 
Kapitanski ]2005], Akhmetiev [2005], and Leal and Pineda ]2008]. 

The principal sources for these formulas are Massey triple products in cohomology, quantum 
field theory in general, and Chern-Simons theory in particular. A common feature of these integral 
formulas is that choices must be made to fix the domain of integration and the value of the 
integrand. 

Our own approach to this problem, initiated in Part 1 and continued here, has been applied 
by Komendarczyk ]2009, 2010] in special cases to derive a higher order helicity for magnetic fields 
whose ordinary helicity is zero, and to obtain from this nonzero lower bounds for the field energy. 
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The key map from ConfsM^ to S^. 

Let X, y and z be three distinct points in Euclidean 3-space . They will typically span a 
triangle there, but are permitted to be colinear, as long as they remain distinct. 

We show a typical configuration below, with the sides of the triangle oriented and labeled as a, 
&, c, with the interior angles labeled as a, /3, 7, and with the orientation of the triangle determining 
a choice of unit normal vector n. 

If the triangle degenerates to a doubly covered line segment, then the sides are still recognizable, 
and likewise the interior angles, with two of them zero and the third 180°. In this case we set n 
equal to the zero vector. 




Figure 1. Triangle and normal vector 



We write [a] = a/\a\ for the unit vector pointing along the oriented side a of our triangle, and 
likewise for the other two sides, and similarly define 

and likewise for the other two pairs of sides. Next, we define a vector in 3-space by the formula 
F{x, y, z) = [a] + [6] + [c] + [b, c] + [c, a] + [a, 6] , 

equivalently, 

( a h c \ 
Fix, y,z) — rr + TTT + TT + (sm a + sm p + sm 7)71 . 
\\a\ \b\ \c\J 

The term [a] = a/\a\ is just the classical Gauss map applied to the vertices y and z, and the 
expression [a] + [b] + [c] is the symmetrization of this. It is a vector tangent to the plane of the 
triangle, and is easily seen to vanish only for equilateral triangles. 

The term [b,c] — b x c/\b\\c\ ~ (sina)n is a dimensionless version of the "directed area" of 
the triangle, and the expression [b, c] + [c, a] + [a, b] is the symmetrization of this. This vector is 
orthogonal to the plane of the triangle, and vanishes only when the triangle degenerates, with the 
vertices lying along a line but remaining distinct. 

It follows that F{x, y, z) is never zero, since it is the sum of two orthogonal vectors that do not 
vanish simultaneously. 

The smoothness of as a function oi x, y and z is apparent from its defining formula, since 
each of its six terms is a smooth function of distinct points. The equivariance of F with respect 
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to orientation preserving isometries of M.^ is similarly apparent, since 

[via)] = f{[a]) and [ip{b), ip{c)] = ip{[b, c]) 

for any rotation ip of M'^, while translations don't change a, b and c. F is likewise insensitive to 
change of scale. 

Let ConfsE^ denote the configuration space of ordered triples of distinct points in M'^. Then 
with the above definition, we have 

F : ConfsM^ ^ - {0} . 
Since the image of F misses the origin of M.^, we may normalize to obtain 

/ = ^ : ConfaM^ ^ . 
Then the map / is also smooth, equivariant as above, and insensitive to change of scale. 
The generalized Gauss map. 

Suppose now that L is a link in M'^ with three parametrized components 

X ^{x{s) : se S^}, Y ^{y{t) : te S^}, Z = {z{u) : u e S^} . 
We define the generalized Gauss map -.T^ ^ S'^ by 

9l{s, t, u) = f{x{s),y{t),z{u)) . 

We regard this map as a natural generalization of the classical Gauss map from the 2-torus to 
the 2-sphere associated with a two-component link in R^. If the link L is smooth, then so is 
the map g^. 

The map is equivariant with respect to the group Isom^ of orientation-preserving isome- 
tries of M"^. That is, if (p is such an isometry, then gip{L) = V ° 9l, where Lp acts on 5^ via its 
"rotational part". In particular, if 1^9 is a translation, then g^{L) — 9l- 

Since the map / is insensitive to change of scale, so is the map g^. 

The homotopy class of g^ is unchanged under reparametrization of L, or more generally under 
any link homotopy of L. The generalized Gauss map is also "sign symmetric" in that it trans- 
forms under any permutation of the components of L by precomposing with the corresponding 
permutation automorphism of multiplied by the sign of the permutation. 

Pictures of the generalized Gauss map. 

In each of the three figures below, we show the vector 

F{x,y,z) = ( + 7TT + A ) + (^™" + sin/^ -Hsin7)n 
\\a\ \b\ \c\J 

attached to a triangle in with vertices at x, y and z. 

Equilateral triangle. 

In this case, (a/|a| -|- b/\b\ + c/\c\) — 0, and so 

F = (sin 60° + sin 60° -I- sin 60° )7i = 3(V3/2)n « 2.6 n, 
as depicted in Figure 2(a). 
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(a) The vector F for an equilateral triangle. (b) The vector F for a 30-60-90 right triangle. 



V 




h 









c 





(c) The vector F for a degenerate triangle. 



Figure 2. The vector F for various triangles 

30-60-90 right triangle. 

In Figure 2(b), we take a to be of length y^, b to be of length 2, and c to be of length 1. The 
vector (a/|a| -f &/|&| -I- c/|c|) is shown above running from a point on the hypotenuse 6 to a point 
on the longer side a, and has length « .52. The quantity 

(sin90° -|-sin60° -f sin30°)n = (1 + \/3/2 + l/2)n w 2.37n, 

and so 

F sa (horizontal vector of length .52) + 2.37 n . 

Degenerate triangle. 

The degenerate triangle shown in Figure 2(c) has a of length 1, 6 of length 2, and c of length 
3. The vector (a/|a| -I- b/\b\ + c/\c\) is then of length 1 lying in the line of the triangle as shown, 
while 

(sin0°+sin0°+sinl80°)n = 0, 
and so F is the unit vector shown. 

Statement of results. 

The first of our two main results gives an explicit correspondence between the Milnor link ho- 
motopy invariants of a three-component ordered, oriented link in and the Pontryagin homotopy 
invariants of its generalized Gauss map. 

Theorem A. Let L be a three- component link in M^. Then the pairwise linking numbers 
p, q and r of L are equal to the degrees of its generalized Gauss map ql ■ — >■ on 
the two-dimensional coordinate subtori of , while twice Milnor's ^-invariant for L 
is equal to Pontryagin's v -invariant for ql modulo 2gcd(p, q, r). 

Each two-dimensional coordinate subtorus of is oriented to have positive intersection with 
the remaining circle factor. 
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We refer the reader to Part 1 for a discussion of Milnor's /i-invariant for a three-component Hnk 
in or , and of Pontryagin's i^-invariant in the special case of a smooth map of a 3-manifold 
to the 2-sphere. We also explained there how to convert Pontryagin's original relative invariant to 
an absolute invariant for maps from the 3-torus to the 2-sphere by comparing with an appropriate 
family of "base maps". 

In Part 1, the long and detailed proof of Theorem A in the spherical case was carried out in 
terms of framed bordism of framed links in the 3-torus. We will capitalize on that effort here by 
building a bridge from the Euclidean to the spherical versions, and cross it to transfer the burden 
of proof from the Euclidean side back to the spherical side, to work already done there. 

Our second main result provides an integral formula for Milnor's /_t-invariant in the special case 
when the pairwise linking numbers p, q and r of L vanish. 

To explain the symbols that appear in that formula, let ui denote the usual area form on S^, 
normalized to have total area 1. Then w pulls back under the generalized Gauss map gL'-T^—^S^ 
to a closed 2-form lul on T^, which we refer to as the characteristic 2- form of L. 

When p, q and r are all zero, it follows from Theorem A that ujl is exact. We then let d~'^{ujL) 
denote any "primitive" of uj^, meaning any 1-form on whose exterior derivative is lul. 

Theorem B. Let L be a three- component link in M.^ whose pairwise linking numbers 
are all zero. Then Milnor's ^-invariant of L is given by the formula 

(1) Mi) = J / d-^{u:L)/\UJL. 

The value of this integral is easily seen to be independent of the choice of primitive c?~^(wl) for 

The geometrically natural choice for d^^{ujL) is the primitive of least L^-norm. It can be 
obtained explicitly by convolving lul with the fundamental solution Lp of the scalar Laplacian on 
, and then taking the exterior co-derivative 5 of the resulting 2-form: 

d^^{u!L) = 5{(p * wl) . 

Details of this construction are presented after the proof of Theorem B. If we make this geomet- 
rically natural choice for cZ~^(wl), then the integrand in the above formula is also geometrically 
natural in the sense that it is unchanged if L is moved by an orientation-preserving isometry of 
M?. 

For comparison with formula (1) above, we recall J.H.C. Whitehead's integral formula for the 
Hopf invariant of a smooth map f : S^, 

Hopf(/)= / d-\f*Lo)M''uj, 
where w is the normalized area form on 5*^, and /*w is its puUback via / to an exact 2-form on 

There are two additional versions of the integral formula for Milnor's /i-invariant given in 
Theorem B, and we present them next. To state these formulas, we again need some definitions. 

Let L be a three-component link in M?, and its characteristic 2-form on T"^. We convert the 
closed 2-form cj^ to a divergence-free vector field Vl on via the usual formula, 

ul{A,B) = {AxB).Vl. 
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for all vector fields A and B on T^. We refer to Vl as the characteristic vector field of L on 
T^. When the pairwise linking numbers p, q and r of L are all zero, the vector field Vl on is 
in the image of curl. 

For the third version of our formula for Milnor's /x- invariant, we need to express the characteristic 
2-form and vector field in terms of Fourier series on the 3-torus. To that end, view ~ S^xS^x 
as the quotient (IR/27rZ)^, and write x — {s,t,u) for a general point there. 

Using the complex form of Fourier series, express 

= ^ (c^ dt Adu + duAds + ds A dt)e'-" ' . 

We compress notation by writing 

dx — (ds, dt, du) , ★dx = {dt A du , du Ads , ds A dt) , 

Using this compression, the formulas for lol and Vl become 

UJL^Y^ c„e'" • ^ ■ * da; and Vl c„e*" " ^ • . 

n n 

Writing = (0,0,0), the coefficient Cq = (cq,Cq,Cq) = since the form ujl is exact, equivalently 
the vector field Vl is in the image of curl. Finally, we express the general Fourier coefficient Cn in 
terms of its real and imaginary parts, 

Cn = + with and bn real. 

Theorem B (continued). Let L be a three- component link in with pairwise linking 
numbers all zero. Then Milnor's ^-invariant of L is also given by the formulas 

(2) ti{L) Vl(x) xVLiy) 'Vy<p{x-y)dxdy 

(3) = Stt^ V a„ X bn . , 

— nM 

where (p is the fundamental solution of the scalar Laplacian on the 3-torus, Vl is the 
characteristic vector field of L, and an and bn are the real and imaginary parts of the 
Fourier coefficients Cn of Vl ■ 

In formula (2) above, the difference a; — y is taken in the abelian group structure on T^, the 
expression Vyip{x — y) indicates the gradient with respect to y while x is held fixed, and dx and 
dy are volume elements on . 

Formula (2) is just the vector field version of formula (1), in which the integral hidden in the 
convolution formula for d~^{uJL) is expressed openly. This formula shows that the Milnor triple 
linking number /i(L) is one-half the helicity of the vector field Vl on the 3-torus . The integrand 
in formula (2) is invariant under the group Isom^ M"^ of orientation-preserving isometrics of M'^. 
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Numerical computation. 

We used Matlab to calculate an approximation to Milnor's /^-invariant, as given by formula (3) 
of Theorem B, for the three-component link L in parametrized by 

x{s) = (2 cos s, 7 sin s, 0) 
y{t) = (0,2cost, 7sint) 
z{u) = (7 sin u, 0, 2 cos u) , 

with s £ [0, 27t], t £ [0, 27r]. u £ [0, 2tt], which is a concrete realization of the Borromean rings with 
/i = — 1 shown in Figure 3. 




Figure 3. Borromean rings 

In particular, we used Matlab to calculate approximations to the Fourier coefficients Cn of its 
characteristic form lul. We used subdivisions of the s, t and u intervals into 256 subintervals to 
approximate the integrals defining the coefficients Cn for —64 < ns,nt,nu < 64. The approximation 
of fj, we obtained in this way was —0.99999997. 

II. Theorem A 

Proof plan for Theorem A. 

We begin with a quick summary of the spherical theory, and its key map fs ■ Gonf^S^ — ?> S"^. 
Next, we discuss inverse stereographic projection h from to S^, and use it to define a map 
H: ConfsM'^ CoufsS*"^. Afterwords, we state and prove the "bridge lemma", which asserts the 
honiotopy commutativity of the diagram 

ConfsS'^ 

H 

ConfsM^ 52 

where the horizontal maps are the "key maps" of the spherical and Euclidean theories. Finally, we 
use the bridge lemma to prove Theorem A. 
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Recollection of the spherical theory. 

The key map ConfaS''^ S"^ in the spherical setting was defined in Part 1 as follows. Let x, y 
and z be three distinct points on the unit 3-sphere in M^. They cannot lie on a straight line in 
M'*, so must span a 2-plane there. Translate this plane to pass through the origin, and then orient 
it so that the vectors a; — z and y~z form a positive basis. The result is an element G{x,y, z) of the 
Grassmann manifold of all oriented 2-planes through the origin in 4-space. This procedure 

defines the Grassmann map 

G: QonUS^ ^ G2M\ 
where ConfaS''^ is the configuration space of ordered triples of distinct points in . 




Figure 4. Map from ConhS-' to G2M^ 

The Grassmann manifold with its natural Riemannian metric is, up to scale, isometric to 

the product S*^ x S"^ of two unit 2-spheres. We will express this by the map tt : — > S*^ x S*^ 
which takes the oriented 2-plane (a, b) with orthonornial basis a, 6 to the point {ba,ab) in x S'^, 
using quaternion notation and conjugation. This gives us two projection maps 7r_|_ and 7r_ from 
G2K^ ^ 5*2, 

7r+ (a, b) = ba and 7r_ (a, b) =ab. 

If the basis a, b is not necessarily orthonormal, then we saw in Part 1 that 

Im(6a) Im(a5) 
TT+{a,b) — J- — and n^{a,b) = 



|Im(6a)| ~' ' ' \lm{ab)\ ' 

We arbitrarily use the first projection 7r+ to define the key map 

fs = TT+oG: ConiaS^ ^ . 
Inverse stereographic projection. 

The corresponding key map in the Euclidean theory is 

\^e\ 

where 

FE{x,y,z) = ( rr + 777 + rr ) + (sin a + sin /? + sin 7)71 
\\a\ \b\ \c\J 

was defined earlier, and where we have added the subscripts to and Fe to signal "Euclidean". 
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The bridge between the two theories will be a map H : ConfaM'^ — >■ ConfsS''^ which makes the 
diagram 



ConfsS'' 

H 



'3 fs 



ConfsM^ — 52 



commutative up to homotopy. We present two versions of H, the first straightforward via inverse 
stereographic projection, and the second homotopic to it but more convenient for our arguments. 

Viewing M'* as the space of quaternions, we regard as the subspacc of purely imaginary 
quaternions, and then use inverse stereographic projection from —1 to provide a diffeomorphism 
h:'M? ^ — {— 1}, which preserves the usual orientations on and S^. 




Figure 5. Inverse stereographic projection h: M? — > — {—1} 
Let q denote a purely imaginary quaternion, thus a point of M?. We compute that 

with the first term on the right being the real part of h{q), and the second term its imaginary 
part. Indeed, a quick check shows that h{q) has norm 1, and that h{q) — (—1) is a real multiple of 
q— (—1), and hence that the points — 1 , and h{q) lie on a straight line. 

The first version of the map H: ConfsR^ ConfaS^ uses inverse stereographic projection on 
each of three points, 

H{x, y, z) = {h{x), h{y), h{z)) . 
The second version, call it H', is defined as follows. Let (x, y, z) be a triple of distinct points in 
R^, and use translation by — z there to move this to the triple {x — z,y — z,0) of distinct points. 
Then apply h to each of the three points in this new triple to obtain 

H'{x,y,z) = {h{x-z),h{y-z),l). 
The maps H and H' : ConfsM"^ Confa^^ are clearly homotopic. 
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Figure 6. The maps H and H' : Confal 
Statement of the bridge lemma. 



Conf,53 



The resuh below will permit us to transfer the burden of proof for our current Euclidean version 
of Theorem A back to its spherical version in Part 1 . 

Bridge Lemma. The map H : ConfaM^ — >■ ConfaS^^ makes the diagram 

ConfaS^ 

I 

H 

ConfaM^ — 



commutative up to homotopy. 

Setup. We avoid the nuisance of normalization by using instead the maps 
Fs : ConfaS'^ ^ - {0} and Fe : ConfaM^ ^ - {0} 

defined by 

Fs{u, V, w) = Im^(ti — 'w){u — wfj 



and 



FE{x,y,z) 



a b c , , . . „ . . 

r-r + TTT + TT + (sm a + sm p + sm 7jn . 
\a\ \b\ \c\ 



At the same time, we replace the map H by the homotopic map H' . So now our job is to show 
homotopy commutativity of the diagram 



ConfaS'^ ^ ^ TO) 3 



H' 



ConfsM^ — — 



- {0} 

- {0} 
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Proof of the Bridge Lemma. 

We start with three distinct points x, y and z in M'^, forming a possibly degenerate triangle with 
sides a — z — y, h — x ^ z and c ^ y — x. Then we begin to compute, 

Fs o H'ix, y, z) = Fs{h{x - z),h{y " 1) 
= Fs{h{b)M~ci)A) 

= Im ({h{-a) - imb) - 1)) . 



We recall the formula 

l-|g|' 2g 
and first substitute —a for g, and then b for q, to get 

/i(-a) - 1 = — 

and 

|2 



;i(6) - 1 = - 



i + w ■ 

Then 

lin[{h{-a) - l)(/i(6)-l)) = CIm((|a|2 + a)(|6p + fe)) 
where the positive real number C is given by 

C = 4(l + |an-i(l + |6|2)-i. 

We keep in mind that a and b, since they lie in M?, are purely imaginary quaternions, and hence 
the sums |ap + a and |6p + b are both quaternions written in terms of their real and imaginary 
parts. 

We recall the following formula about quaternion multiplication, 

Im(gig2) = Re(gi) Im(g2) + Im(qi) Re(g2) + Im((7i) x Im((72) , 
using the vector cross product in the 3-space of purely imaginary quaternions. 
Applying this formula, we get 

Im((|a|2 + a)(|6p + fe)) = \a\% + a\b\'^ + a x b . 
Then, stringing together the above computations, we have shown that 

Fs o H'{x, y, z) = C{\a\% + a\b\^ + a x b) 

with 

C = 4(l + |an-i(l + |6|2)-i. 

Since we are focusing on homotopy of maps into M.^ — {0}, the strictly positive quantity C is 
irrelevant, and we hide it from view by recording the homotopy 

(1) Fs o H\x, y, z) ~ |ap6 + a\b\^ + axb. 

It remains to show that this expression is homotopic in M'^ — {0} to 

(2) FEix,y,z) = ( + TTT + A ) + (sina + sin/3 + sin7)n. 

\\a\ \b\ \c\J 

The right side of (1) can be rewritten as 
(1') {\a\^b + a\b\^)+{\a\\b\sm^)n 
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for convenience of comparison with (2). In each case we have the sum of a vector parallel to the 
plane of the triangle xyz and a vector orthogonal to it. 

As long as the triangle is non-degenerate, the components in (!') and (2) orthogonal to its plane 
are both strictly positive multiples of the unit normal vector n. Hence (!') and (2) are vectors 
which both lie in the same open half-space of 'M? , and so the line segment between them misses 
the origin. Thus Fs o H' is homotopic to Fe in such cases. 

So the issue now is, what happens when the triangle degenerates? In such a case, (!') and (2) 
reduce to their tangential components, 

(1') |a|2&-fa|&|2 
and 

a h c 

H + H^R- 

Both of these vectors point along the line of the degenerate triangle, and we must check that they 
always point the same way, so that the line segment between them again misses the origin. 



X c ^ y a ^ z 




Figure 7. Degenerate triangles 



There are three cases, according as which vertex is between the other two. They are shown in 
Figure 7, in which the line of the degenerate triangle is turned so that the two shorter sides point 
to the right. 

In each case, the vector a/\a\ + b/\b\ + c/\c\ is a unit vector pointing to the right. Furthermore, 
in all three cases, the vector |ap6-|-a|&p points in the same direction as its rescaling -|-a/|ap, 
which is a nonzero vector pointing in the same direction as the shorter of the two vectors a and b, 
and this is also to the right. It follows that in every case, non-degenerate or degenerate, the line 
segment between the vectors (!') and (2) misses the origin. 

Hence the maps Fs o H' and Fe from ConfaM^ — ?> M"^ — {0} are homotopic, completing the proof 
of the Bridge Lemma. 

Proof of Theorem A. 

Let L be a three-component link in M'^, let h: M.^ 5"^ — { — 1} be inverse-stereographic projec- 
tion, and let h{L) be the resulting three-component link in . Since h is an orientation-preserving 
diffeomorphism, the Milnor invariants p, q, r and /i for L match those for h{L). 

The Euclidean generalized Gauss map : ^ for L is given by 

9l{s, t, u) = fs {x{s), y{t),z{u)) , 
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while the spherical generalized Gauss map gh{L) : ^ for h{L) is given by 

9h{L){s,t,u) = fs{h{x{s)),h{y{t)),h{z{u))) 
= fs ° H {x{s),y{t),z{u)) . 
According to the Bridge Lemma, the maps fs and fs ° H : Conf sM'^ — !• are homotopic. 

It follows that the maps and gh{L)'- are also homotopic. Hence the Pontryagin 

invariants p, q, r and v for g^ match those for gh(L)- 

Then the correspondence between the Milnor invariants for L and the Pontryagin invariants for 
gL, as asserted in our current Euclidean version of Theorem A, follows from the correspondence 
between these invariants for h{L) and gh(L), as asserted in the spherical version of Theorem A, 
which was proved in Part 1. 

This completes the proof of the Euclidean version of Theorem A. 

III. Theorem B 

Proof plan for Theorem B. 

Let L be a three-component link in Euclidean space M.^ , with pairwise linking numbers all zero. 
Theorem B gives three explicit formulas for the triple linking number (Milnor ^-invariant) of L: 

(1) ^l{L) = ^J^J-\oJL)AoJL 

(2) VL{x)xVL{y) •^yV{x~y)dxdy 

(3) = ^ a„ X b„ . , 

n^O ' ' 

using the notation defined after the two statements of Theorem B. 

We saw in Part 1 that the spherical version of Theorem A implies the spherical version of 
Theorem B, and we show below that the same implication holds for the Euclidean versions here. 
We will give details only for the proof of formula (1), and refer the reader to Part 1 for the 
derivation of formulas (2) and (3). 

The first step will be to give an explicit formula for the characteristic 2- form W]^ of the link L. 

As mentioned earlier, the geometrically natural choice for d~^{ujL) is the primitive of least 
L^-norm, which can be obtained explicitly by convolving lul with the fundamental solution ip of 
the scalar Laplacian on T'^, and then taking the exterior co-derivative S of the resulting 2-form: 
d^'^{ujL) = 6{ip * ujl). 

So we will give the expression for this fundamental solution ip, and then show how formula (1) 
follows from J.H.C. Whitehead's integral formula for the Hopf invariant. 

An explicit formula for the characteristic 2-form ljl. 

We start with a three-component link L in M.^ with components 

X ^{x{s) : se S^}, Y ^{y{t) : te S^}, Z = {z{u) : u G S^} , 
and recall the figure 



15 



DeTurck, Gluck, Komendarczyk, Melvin, Nuchi, Shonkwiler and Vela-Vick 




with 

a{t, u) — z{u) — y{t), b{s, u) = x{s) — z{u), c(s, t) — yit) — x{s) , 
and the formula 

F{x, y, z) = [a] + [b] + [c] + [h, c] + [c, a] + [a, b] . 
With mild abuse of notation, we write 

F{s,t,u) ^ F{x{s),y{t),z{u)) . 
The generalized Gauss map ql ^ of the link L is then given by 

F(s,t,u) 



9L{s,t,u) 



\F{s,t,u)\ 



Let uj be the Euclidean area 2-form on the unit 2-sphere S"^ C M"^, normalized so that the total 
area is 1 instead of An. If P is a point of S'^, and A and B are tangent vectors to S'^ at P, then 

ojp(A,B) = —{A X B) • P. 
All 

This 2-form uj on S"^ extends to a closed 2-form TU on M?" — {0} given by 

, , , (AxB) • P 

which is the puUback of cj from S*^ to - {0} via the map P ^ P/\P\. 

Hence the puUback g*^(jj of w from S"^ to via — Fl\F\ is the same as the puUback F*u: of 
id from M? ~ {0} to via F. Write 

g*j^Lj = F*Ld = p{s, t, u) dt A du + q{s, t, u) du A ds + r(s, t, u) ds A dt . 

Then we have 

p{s, t, u) = F*uj{dt, du) = uj{F^dt,F^du) 
{F,xFu)-F 

where the subscripts on F denote partial derivatives, and likewise for q(s,t,u) and r{s,t,u). 
Therefore, the characteristic 2-form of the link L is 

= — ^ ((Ft X Fu- F)dtAdu+ (F„ x F, • F) du A ds + (F, x Ft • F) ds A dt) . 
47r F P V / 
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Proof of Theorem B, formula (1). 

Let L be a three-component link in Euclidean 3-space M.^ with pairwise linking numbers p, q 
and r all zero. 

By the first part of Theorem A these numbers are the degrees of the Gauss map : T"^ — > S*^ 
on the two-dimensional coordinate subtori. Since these degrees are all zero, is homotopic to a 
map g: ^ which collapses the 2-skeleton of to a point: 

where a is the collapsing map. 

By the second part of Theorem A, Milnor's /z-invariant of L is equal to half of Pontryagin's 
z^-invariant of g^, which in turn is just the Hopf invariant of /: 5"^ — > S^, 

^Ji{L) = \y{gL) = \^ov^{.f). 
We can thus use J.H.C. Whitehead's integral formula for the Hopf invariant, as follows. 

Let be the area 2-form on S"^ , normalized so that ^g^^ = 1- Its puUback /*a; is a closed 
2-form on , which is exact because iJ^(5''^;E) = 0. Let d~'^{f*uj) indicate any smooth 1-form 
on whose exterior derivative is f*uj. Then, as noted earlier. Whitehead showed that the Hopf 
invariant of f is given by the formula 

Hopf(/)= / d-\rLo)Ar^, 

the value of the integral being independent of the choice of the 1-form dr^{f*uj). 

Pulling the integral back to via the collapsing map a: ^ yields the formula 

thanks to the fact that (^l is homotopic to g — / o cr, and recalling that uj^ — <?2'^- Since 
^(L) = \v{gL), we get 

completing the proof of formula (1) of Theorem B. 

The geometrically natural choice for d~^(a;£,). 

After stating Theorem B, we indicated that the geometrically natural choice for d^^{ujL) is the 
primitive of least i^-norm. It can be obtained explicitly by convolving uj^ with the fundamental 
solution ip of the scalar Laplacian on T^, and then taking the exterior co-derivative S of the resulting 
2-form: 

d^^{ujL) = S{ip * wl) . 
If we make this geometrically natural choice for d^^(cji), then the integrand 

d^^{ujL) A iOL 

in the above formula for the triple linking number ii{L) is also geometrically natural in the sense 
that it is unchanged if L is moved by an orientation-preserving isometry of M.^. 

We give a hint of the details, extracted from Part 1, and refer the reader there for proofs of 
Propositions A and B below. 
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Proposition A. The fundamental solution of the scalar Laplacian on the 3-torus 
— (]R/27rZ)^ is given by the formula 



The function ip is C°° at all points x 6 T^ except 0, where it becomes infinite. 

In the above formula, n denotes a triple of integers. 

Although this formula for if is expressed in terms of complex exponentials, the value of ip is real 
for real values of x because of the symmetry of the coefficients. Figure 8 shows the graph of the 
corresponding fundamental solution 



of the scalar Laplacian on the 2-torus T'^ ~ x , summed for |n| < 15, and displayed over the 
range [— Stt, Stt] x [— Stt, Stt]. 



The graph in Figure 8 looks like a "Morse function" with infinite maxes at the lattice points, 
saddles in the middle of the "edges", and mins at the center of the fundamental domains. 

Proposition B. // w is any exact differential form on with C°° coefficients, then 



is a C°° differential form satisfying da. — u). Furthemore, if da — U3 as well, then 
\a\L2 < |a|i,2, with equality if and only if a — a. 






Figure 8. Fundamental solution of the Laplacian 



a — d{ip * Lj) 



18 



Generalized Gauss maps and integrals for three-component links II 



Epilogue 

Where do the generalized Gauss maps come from? 

In the spherical theory, the generalized Gauss map ql'-T^ ^ S"^ comes from the key map 

fs : Conf35=^ Grassmann map G ^^^4 , ^2 

via the substitution 

fi(z,(s, t, u) = fs {x{s),y{t), z{u)) , 

while in the Euclidean theory it comes in the same way from the unit normalization /e of the key 
map 

Fsix, y, z) = [a] + [h] + [c] + [&, c] + [c, a] + [a, b] . 
But where do these key maps come from? 

In the spherical theory, we saw in Part 1 that the configuration space ConfaS^ deformation 

retracts to a subspace difTeoniorphic to x S"^, and the key map fs there is an S'(9(4)-equivariaiit 
version of this deformation retraction, followed by projection to the S*^ factor. 

In the Euclidean theory, we face two complicating features: the configuration space ConfsM^ 
is more challenging - it deformation retracts to a subspace diffeomorphic to a nontrivial S"^ V S"^ 
bundle over 5^ - and the group Isom"*" of orientation-preserving isometries of is non-compact. 

If we were not seeking a generalized Gauss map which is geometrically natural in the sense of 
being Isom"'" M^-equivariant, we could simply define the key map fs to be the composition 

ConfaK^ ^L_. Confs^^ — 5^ , 

where the "inclusion" H : ConfsM"^ — ConfaS''' was defined earlier via inverse stereographic pro- 
jection. This definition of fs is far from being Isom"*" M^-equivariant, and the resulting generalized 
Gauss map would suffer from the same defect, and so lose its applicability to problems in fluid 
dynamics and plasma physics. 

What we did instead was to consider the map H' : ConfaM'^ — > Confs^^ which first took a 
triple (.T, z) of distinct points in M"^ via translation to the triple (x — z,y — z,0) "based" at the 
origin in M.^, and then via inverse stereographic projection h to the triple {h{x — z), h{y — z), 1) of 
distinct points based at the identity in S^. 

Inverse stereographic projection for such based triples is 50(3)-equivariant, and leads to the 
map Fs o H': ConfsM^ — > - {0} which up to scale takes 

(a;, y, z) i — «- |a|^6 + a|6|^ + a x 6 , 

thanks to our earlier computation, where a = z — y and h = x — z. This map is Isom"*" M^- 
equivariant, and so is the resulting projection to S"^, but at the cost of losing scale-invariance and 
"sign symmetry" in the three points x, y and z. 

A little artful play led to the alternative formula 

Fe{x, y, z) = [a] + [b] + [c] + [b, c] + [c, a] + [a, b] , 

which is still Isom"*" ]R"^-equi variant, but now also scale-invariant and sign symmetric, and at the 
same time, thanks to the Bridge Lemma, homotopic to Fs o H'. 

This is the origin of the key map Fe : ConfsR^ ^ — {0} and the resulting Euclidean version 
of the generalized Gauss map g^: ^ S'^. 
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